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Abstract. - When a corrosive solution reaches the limits of a solid sample, a chemical 
fracture occurs. An analytical theory for the probability of this chemical fracture is proposed 
and confirmed by extensive numerical experiments on a two dimensional model. This theory 
follows from the general probability theory of extreme events given by Gumbel. The analytic 
law differs from the Weibull law commonly used to describe mechanical failures for brittle 
materials. However a three parameters fit with the Weibull law gives good results, confirming 
the empirical value of this kind of analysis. 



Chemical etching of disordered solids is an important technological problem that 
presents as well interesting questions in the theory of random systems Strong etching 
solutions will eventually lead to the fracture of a finite sample, an event defined as "chemical 
fracture" . In this letter we present a theory for the statistical behavior of this specific fracture 
mechanism. We show || that the deep nature of the chemical fracture statistics relies on 
the probabilistic theory of extreme events due to Gumbel [Q. Interestingly it is found that 
our chemical fracture statistics are not practically distinguishable from the Weibull statistics 
empirically introduced to fit mechanical fracture statistics [El . 

The problem is studied using a simple two-dimensional corrosion picture || inspired by an 
experimental study of pit corrosion of aluminum films . The model describes the chemical 
etching of a random solid by a finite volume of solution. It predicts that the etching process 
stops spontaneously on a fractal liquid-solid interface as observed experimentally. 

The two dimensional system is described in Fig. ED: The solid is made of lattice sites 
exhibiting random "resistances to corrosion" €E [0, 1] uniformly distributed. It has a width 
L and a depth Y. At any time t the "etching power" of the solution is proportional to the 
etchant concentration C(t) : p(t) — TC(t). The solution has a volume V and contains an 
initial number N et (0) of etchant molecules. Therefore p(0) = C(0) = N et (0)/V using T = 1. 
Hereafter we choose p(0) > p c , where p c is the percolation threshold of the lattice. 

The solution is initially in contact with the solid through the bottom boundary y = 0. 
At each time-step t, all surface sites with T{ < p(t) are dissolved and a particle of etchant is 
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Fig. 1 - Sketch of the etching dynamics in a square lattice: the sites 2, 3, 5 are etched at the first 
time-step as their resistances are lower than p(0). At the same time the number of etchant particles 
in the solution decreases by 3 units, and a new part of the solid is uncovered. The process is then 
iterated. 



consumed for each corroded site. Hence, the solution concentration progressively decreases. 
For an arbitrarily large sample, the corrosion process always stops spontaneously at a certain 
time tf when the etchant concentration p(tf) is still finite [||, reaching a maximal depth 
However, when Y is finite, the corrosion can reach the bottom of the sample before tf, 
fracturing the solid into two disconnected parts. This "chemical fracture" is a stochastic event 
which depends on the realization of the random resistances r*. 

Hereafter we call V(V) the fraction of samples broken by a volume V of solution, keeping 
the other system parameters L, p and Y fixed. Its measure is shown in Fig. |[ The solution 
volume represents the applied chemical force for given p . Obviously the fracture event is 
directly related to the value of the maximal depth yM relative to the sample depth Y: if 
Vm < Y then the fracture does not occur, otherwise the sample will be broken in two distinct 
pieces. Hence, in order to compute the fracture statistics we have to study the probability 
distribution of yu in a sample with an infinite depth. 

Before proceeding, it is important to remark the differences between the chemical fracture 
statistics and the distribution of the maximal corrosion depth. The chemical fracture statistics 
V(y) is a parametric curve which gives the fraction of similar solid samples fractured for a 
given value of V. Hence V(V) is not a distribution function in the sense of probability theory. 
In particular V is not a random variable but an external parameter. 

On the contrary, the distribution of the maximal corrosion depth is the probability distri- 
bution function (DF) of the random variable yu, representing the maximal corrosion depth 
for an (infinitely deep) sample. This is a genuine probability distribution function and, accord- 
ingly, its derivative is the probability density function (pdf) of the variable xjm- Obviously 
the two quantities are related, but of a different nature. In the following we will show how it 
is possible to recover the first, via a probabilistic theory of the second. 

The fracture statistics V(V) of Fig. 2 has been measured numerically taking a sample of a 
fixed depth and performing several corrosion processes for each value of the solution volume 
V, and recording the fraction of broken samples as a function of V (see Fig. |^). 

On the other hand, the probability distribution of the maximal corrosion depth have 
been measured fixing V and recording the maximal corrosion depth during several corrosion 
processes on large systems. A different choice of V give a different distribution, but a universal 
scaling form can be recognized. As shown in Fig. || the numerical DF P(zm < u) for the 
reduced variable zm = (vm — (um))/^ (where (j/m) is the average value of yu and S its 
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Fig. 2 - Chemical fracture probability P (vm > Y) as a function of the volume V which measures the 
"chemical force". Small circles are the numerical estimated probabilities performing 1000 runs for 
each value of V applied on a solid of sizes L — 1000 and Y = 500 (po — 0.7). Diamonds represent 
the fit of the fracture data with a Weibull law of parameters Vo = 1.15 • 10 5 , Vi = 4.08 • 10 , and 
m = 3.89. The line is the result of the fit using Eq. (lid), as explained in the text (the fitted parameter 
A = 0.3875). 



standard deviation) is the same for several choices of V. 

Since ijm is an extremal random variable, we compare this numerical distribution with the 
standard (i.e. with zero mean and unitary variance) Gumbel distribution function H(u) [Q 
for extreme events: 

H(u) = e-e-^ ? (1) 

where a ~ -0.5772 and b ~ Vl.64493. 

We recall that the Gumbel distribution is the canonical distribution for the maximal value 
of a set of random variables in the same manner that the Gaussian distribution is the canonical 
distribution for the sum. More precisely if {xi, xn} is a set of independent stochastic 
variables with identical distribution decaying faster than any power law for large values, then, 
in the large N limit, the maximum value among them can be shown to be Gumbel distributed. 
The collapse of our data on the Gumbel distribution is remarkable, first because there are no 
adjustable parameters, and second because the points on the interface are correlated to some 
extent. 

In the following we discuss why the physics of the etching model determines such a Gumbel 
behavior. As a second step we will derive from it the theoretical form of V(V). This last task is 
accomplished considering the probability P(hm > Y) and introducing the explicit dependence 
of (dm) and S on V (with po, Y and L fixed), i.e.: 



V(V) = 1 — exp< — exp 



\ (Vm){V)-Y 

— a 



(2) 



Such explicit dependence is obtained using the known scaling results of this particular etching 
model. 

We now proceed with the derivation of the Gumbel law for this etching model. First 
we have to explain why the corrosion front can be reduced to a set of independent random 
variables. At the beginning of the dynamics ||, p(t) follows the exponential law p(t) = 
Po exp(— t/r) with r = V/L. The corrosion front advances layer by layer up to approximatively 
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Fig. 3 - The numerical integrated probability distribution for reduced ym (zero mean and unitary 
variance) is compared with standard Gaussian and Gumbel distributions. In the main figure, the 
simulation results collapse on the Gumbel distribution given by Eq. (Q) (hidden by data points). In 
the insets we show the squared difference between data and the Gaussian and Gumbel distributions. 
The data refer to simulations performed for L = 3000, 5000 and iV et (0) = 5 ■ 10 6 , 1 ■ 10 7 , 2 • 10 7 , 5 • 10 7 
(1000 different dynamics for each choice of L and V, po = 0.7). 



t = t c when p{t) — p c and reaches a depth yu n . ~ V/L. After this period the corrosion front 
becomes very irregular and finally stops at t = tf. At tf the etching power Pf = p(tf) is 
slightly smaller than p c and the final corrosion front is fractal with dimension Df = 7/4 up 
to a characteristic width a < L. This situation is displayed in Fig. ||. As shown in ||, this 
phenomenon obeys the scaling laws of Gradient Percolation || where the role of the gradient 
is played by the ratio L/V. This implies that 

• a - (L/V)- 1 /^ 

• a can be seen as a percolation correlation length. 

The total final corrosion front, shown in Fig. ||, can then be considered as a juxtaposition of 
N ~ L/a nearly independent fractal boxes of lateral width (for a similar recognition of 



effective independent variables in extreme statistics see 10 ). Inside each box k, there is a 



, . (k) 

point of maximal penetration of the front y M . Since such points belong to different boxes, 

(k) 

the values y y M ' are, by construction, identically distributed and independent random variables. 
The extreme position of the front is the maximal value between such a collection of random 
variables. This is why a Gumbel distribution is observed. 

(k) 

More precisely due to the underlying percolation phenomena, the DF for each y M has 
exponential tail with the same characteristic scale ^a: 

P (vm >v)~ e- V -^ for y - (y$) > a, (3) 

where c is a constant of order 1 (see also some recent results on the size of sub-critical 
clusters [0). In this case, the theory of extreme statistics 0] imposes that, for N ^> 1, yM is 
asymptotically Gumbel distributed, i.e. P [(j/Af — (um))/^) < u] = -ff(w) and that: 

(Vm) = (yj^ + ElogJV (4) 



£= lea. (5) 
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Fig. 4 - Final corrosion front: It is composed of TV ~ L/er independent regions of size a (a can be 
considered as the correlation length of the system). The extremal front position yu is indicated. It 
is the max between the N independent values of y M of each region. 



It is then necessary to test numerically Eqs. (|J) and (|j5|) (and specifically the logarithmic 
dependence on N, i.e. on L/a), in order to confirm this arguments. In the situation described 

in Fig. |i| one observes that (yff^ is given by the average depth yj of the final corrosion 

front (average taken over all the final front sites) plus a positive shift of order a. In turn, the 
average depth yj (a directly measurable quantity in simulations) is given by the depth yu n . 
reached during the linear part of the corrosion in addition to a shift (sub-leading) again of 
order a. Therefore from Eq. M) we can write: 



(Vm) - Vf 
Uf ~ Vlin. 
yiin. 

s 



a[l + ki In (L/cr)] 

17 

V/L 



(6) 
(7) 
(8) 
(9) 



where k± is a coefficient of order one. Recalling that a ~ {V/L)^ 7 , we can substitute (j/m) 
and E as functions of V in Eqs. (JsJ) and (^) in order to obtain the real functional dependence 
of V(V) on V apart from numerical scaling coefficients. 

The scaling behaviors are themselves confirmed by the extensive numerical simulations 
shown in Fig. Note that if L/a is kept constant (i.e. if L ~ I/ 4 / 11 ) then [(dm) — yj\ ~ 
{L/V)^ 1 / 01 . The measured exponent 0.59 ± 0.02 is consistent with this prediction. In the 
same figure, we also report the scaling behavior of the standard deviation E and of y/. These 
fits confirm that E ~ a and that yf can be written as the sum of yu n . and a shift of order a. 
Finally, the linear logarithmic dependence on L for fixed a (i.e. for fixed L/V), which is the 
essential confirmation of Eq. (^), is shown in the inset. It has also been numerically confirmed 
that E and yf does not depend on L for fixed sigma (not shown in the figure). 

To summarize, Eq. (0) and Eqs. (||- 1|) give the theoretical form for V(V), once the propor- 
tionality coefficients in Eqs. (|6|-0) and the coefficient k\ have been determined. For simplicity, 
it is convenient to rewrite Eq. (Q) as follows: 



b(y f - Y) 



(10) 



The dependence on k±, described in Eq. (^|), is now contained in the coefficient A. As a closure 
test of our analysis we choose to use this form to fit against A the numerical fracture statistics 
(i.e. the fraction of simulation runs breaking the sample as a function of V) [Q. As shown 
in Fig. 0, the theoretical form of V(V) and the simulation are in excellent agreement. 
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Fig. 5 - Scaling behavior of yf, [{um} — Vf] an d S. The data in the main graph refer to simulations 
performed with a constant L/a (range of parameters: from L = 522 and V = 5 ■ 10 8 to L — 11550 
and V = 5 • 10 s ). The numerical fits are in accordance with the theoretical behaviors. The inset 
shows the logarithmic dependence of (j/m) on L for a given a (parameters range from L = 1205 and 
V = 120500 to L = 80000 and V = 8 ■ 10 6 , chosen to have a constant a). In both cases the data 
correspond to 1000 different dynamics for each choice of L and V, po = 0.7. 

In summary, we have studied the chemical fracture due to the extremal propagation of a 
corrosion front in a two-dimensional scalar etching model. It has been shown both theoretically 
and numerically that the statistics of the maximal depth reached by the solution is given by 
the well known Gumbel extremal distribution function. From this it has been possible to 
extract theoretically the probability of a chemical fracture. 

Before concluding, we believe that it could be interesting to compare our results with the 
commonly analysis used for mechanical fracture statistics. In this case, even though a widely 
accepted theoretical frame is still lacking, the standard empirical fit is often performed using 
a Weibull law 0. In our chemical fracture model, this comes to try to fit the datas with: 



where Vq (the minimal stress to have a finite fraction of fractures), V\ and m are suitable 
parameters. 

The result, shown in Fig. 2, is that the chemical fracture statistics can be fitted very 
accurately with a Weibull law. We have also verified (not shown here) that the Weibull low 
probability tails (the test generally used by engineers) fits nicely with our data. On one side, 
we believe that this is astonishing, since, despite the simple scalar nature of the model, we 
recover, for our purely chemical process, the same empirical statistics found for mechanical 
failures of brittle materials, generally observed in more complex (vectorial or tensorial) frames. 
On the other side, it is important to remark that in the present case the Weibull law is only 
an empirical fit, because the real analytical form of the chemical fracture statistics is given 
by Eqj^. Moreover a Weibull law, which is usually related to the presence of power law 
distributed flaws, has no reason to emerge here, as Eq.|^ results from a theoretical calculation 
of the model and directly connected to known results in percolation theory. 

As a final remark, one could suggest that the correspondence between the Weibull fit 
and the underlying Gumbel distribution (first suggested by Q) found here is linked more 




with V > Vq , 



(11) 
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generally to Gradient Percolation. Gradient Percolation situations can also be present in 
solids which have never been etched ]l4| . The essential ingredients of our theory are (i) the 
percolation aspect linked to randomness of the system, and (ii) the existence of a gradient. In 
real systems gradients may exist in the sample, for example the ion concentration gradient in 
glass fibers used in communication optics. But it may also appear in experiments, for instance 
in experimental studies of ceramic fracture. In classical flexure experiments the distribution of 
stress is non uniform [ |D| . It remains to be determined specifically if such gradient mechanisms 
are responsible for some of the ubiquitous apparent Weibull statistics. 

Our two dimensional results could be generalized to higher dimensions, according to |l6| ] . 
Therein, field theoretic arguments show that the corrosion dynamics towards the final state 
can be described as a self-organized absorbing phase transition towards the critical phase of 
percolation in any dimension. 

* * * 

We are grateful to M. Dejmek and S. Zapperi for useful comments and we acknowledge 
the support of the European Community TMR Network ERBFMRXCT980183. 
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